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INTRODUCTION 


This  paper  explores  the  mathematical  relationships  between  beam  vibration  and  imposed 
launch  loading.  Launch  loads  include  forces  from  down-barrel  coil  reaction,  traditional  gas 
-  propellent  pressures,  and  projectile-barrel  interaction  loads.  Ben6t  Laboratories'  proven  Uniform 
Segments  Method  (USM)  gun  vibration  model  is  used  to  analyze  and  simulate  firing  dynamic 
response  of  a  60-mm  hybrid  cannon.  The  cannon  was  designed  and  constructed  at  Benet 
Laboratories  in  collaboration  with  an  electrical  engineering  group  from  Polytechnic  University  at 
Brooklyn,  NY. 

BACKGROUND 

Great  strides  have  been  made  in  recent  years  to  identify  the  dominant  loads  that  cause 
beam-type  vibrations  in  cannons  during  ballistic  operation.  These  motions  infringe  upon  shot 
accuracy,  because  at  projectile  disengagement,  the  muzzle's  kinematic  state  may  compromise  the 
projectile's  intended  flight  path.  In  the  case  of  gas  guns,  the  dominant  loads  that  cause  these 
motions  have  been  identified  as  interactions  with  the  cannon's  supports  and  eccentrically-applied 
inertia  loads  due  to  non-centered  masses  attached  to  the  barrel  (breech,  bore  evacuator,  muzzle 
brake)  or  non-symmetric  loads  due  to  take-up  of  slack  in  support  and  recoil  components. 
Significant  muzzle  motions  have  been  modeled  and  reproduced  in  tests  when  these 
conditions  exist. 

The  incorporation  of  electromagnetic  propulsion  at  the  muzzle  end  of  a  traditional  gas 
gun  invokes  a  new  load  on  the  structure.  This  loading  is  manifest  as  a  coaxial  reaction  force 
imposed  on  the  muzzle  end  of  the  barrel.  It  may  be  modeled  as  ^  follower  force  that  moves  with 
the  barrel  during  vibrational  deformation,  such  that  the  reaction  load  is  always  tangent  to  the 
deformed  center  line  of  the  barrel  at  the  location  of  the  electromagnetic  accelerator.  This  force 
places  the  barrel  in  a  state  of  columnar  compression  that  tends  to  exacerbate  the  transverse 
vibrations  of  the  Gun  system  during  launch.  This  in  turn  may  lead  to  an  undesirable  increase  in 
round-to-round  dispersion. 

This  paper  presents  a  comprehensive  model  of  hybrid  gun  dynamics  that  allows  computer 
simulation  of  various  parametric  configurations.  This  will  enable  the  engineer  to  include  design- 
for-accuracy  considerations  and  provide  the  opportunity  to  modify  barrel  stiffness,  axial  position 
of  the  electromagnetic  accelerator,  and  mount  geometry  to  avoid  shot  ejection  at  a  time  of 
significant  structural  deflection. 

ANALYSIS  AND  MODELING  OF  PROFILE-DEPENDENT  LOADS 

In  this  paper  the  gun  vibrations  are  analyzed  from  the  perspective  of  transverse  beam 
dynamics.  To  achieve  this,  the  gross  rearward  motion  of  the  gun  system — due  to  the  rigid-body 
recoil  of  the  barrel — must  be  decoupled  from  its  transverse  vibrations.  In  this  section,  equivalent 
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transverse  forcing  functions  will  be  developed,  and  a  solution  provided  based  on  modal 
superposition. 

Recoil  Inertia  Load 

Recoil  is  the  reaction  force  placed  on  the  barrel,  by  Newton's  third  law,  that  balances  the 
launch  loading  imposed  on  the  projectile.  The  recoil  direction  of  a  gun  is  a  function  of  the  mount 
within  which  it  resides.  Current  design  philosophy  specifies  that  rail  and  channel  sliding  surfaces 
be  used  in  future  gas  gun  designs.  This  method  of  support  attenuates  gun  curvature  within  the 
mount,  maintaining  a  straight  pull  that  tends  to  “snap”  the  overhanging  muzzle  end  of  the  gun. 

An  inertia  couple  occurs  that  is  proportional  to  the  slope  of  the  deflection  curve  and  recoil 
acceleration.  This  load  is  expressible  for  a  uniform  beam  segment  as  follows: 

fr(t^x,yjy'^)  =  pa[{L-x)y']'  .jx 


where: 


fr 

a 

t 

L 

X 

y 

Vd 

ys 


p 
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recoil  inertia  load  per  unit  length 
a(t)  =  recoil  acceleration 

time 

length  of  gun 

spatial  coordinate  along  bore  axis 

=  transverse  displacement  of  bore 
yj(x,t)  =  dynamic  displacement  of  bore 

y^Cx)  =  static  displacement  of  bore  (caused  by  gravity 

droop,  manufacture,  etc.) 
p(x)  =  linear  density  of  gun 

d/d\ 


Pressure  Curvature  Load 


Due  to  the  nature  of  curvature  within  beam-type  structures,  diametrically  opposite  bore 
surfaces  possess  differing  areas.  A  pressure  load  acting  within  the  bore  tends  to  straighten  the 
gun  tube  because  the  area  of  the  concave  surface  is  greater  than  its  convex  counterpart.  This 
phenomenon  is  often  referred  to  as  the  “Bourdon  effect.”  This  type  of  load  travels  behind  the 
projectile  and  is  proportional  to  the  second  derivative  of  the  deflection  curve.  The  expression  for 
this  load  function  is 
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where: 


pressure  curvature  load  per  unit  length 
bore  area  of  tube 

=  propellant  gas  pressure 

Heaviside  unit  step  function 

H  =  1  for  0  <  X  i  Xp 

H  =  0  for  X  >  Xp 

Xp(t)  =  axial  location  of  projectile 


The  accelerating  projectile,  although  considerably  less  massive  than  the  tube,  can  exert  a 
significant  transverse  force  when  it  is  constrained  to  travel  along  a  curved  path.  This  load,  which 
contains  transverse,  Coriolis,  and  centrifugal  accelerations,  is  modeled  as  a  point  load  applied  at 
the  instantaneous  axial  location  of  the  projectile  and  is  expressed  as  follows:^"’ 

=  (3) 

-%[y  +  g  +  2ipy  4  Xp2y // ]  6(xp  -  X ) 


fc 

= 

Pb 

H(Xp-x) 


Projectile  Trajectory  Load 


where: 


-X) 


Electromagnetic  Load 


projectile  trajectory  load  per  unit  length 
acceleration  of  gravity 
projectile  mass 

Dirac  6  function,  in  units  of  length'^ 
6(Xp  -  x)  =  0  for  X  Xp 

J6(x^-x)dx  =  1 
d/dt 


The  reactive  loading  placed  on  the  barrel  by  the  electromagnetic  augmentation  of  a  short 
stack  of  coils  may  be  modeled  as  &  follower  force  positioned  at  the  geometric  center  of  the  coil 
stack.  As  such,  this  load,  which  is  computed  off-line  using  a  finite  element  method,  is  oriented 
along  the  axis  of  the  coil  pack  and  is  subjected  to  the  lateral  and  angular  deflection  of  the  barrel 
at  that  point.  This  force  may  be  resolved  into  two  orthogonal  components:  parallel  to  the 
undeformed  center  line  axis,  and  perpendicular  to  the  undeformed  axis.  The  perpendicular 
component  may  be  directly  added  to  the  lateral  forcing  function  at  the  application  point.  The 
parallel  component  pushes  the  entire  gun  system  forwards  or  rearwards  (a  decoupled  rigid-body 
loading)  and  effects  a  distributed  moment  along  the  length  of  the  barrel  (a  distributed  transverse 
loading)  expressible  as:  M(t,x)  =  feJt)cos(-/(x)lj^.^^J[y(x)  -  y(x)lj^-xeJ-  In  this  expression,  the 
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force  and  cosine  terms  resolve  the  parallel  component  of  the  time-dependent  point 
electromagnetic  loading.  The  term  in  the  brackets  represents  the  effective  lateral  moment  arm 
between  any  point  along  the  barrel  and  the  point  of  electromagnetic  load  application  at  In 
addition  to  the  distributed  moment,  the  Euler  beam  reaction  to  the  parallel  component  effects  a 
lateral  point  load  contribution  that  cancels  the  perpendicular  component  under  the  small  angle 
approximation.  (For  elaboration  of  follower  forces,  see  Reference  [5]).  Explicit  derivation  of 
this  load  function  is  explained  below.  The  equivalent  lateral  load  of  a  distributed  moment  of  the 
Euler  beam  is  its  second  spatial  derivative.^*^  The  resulting  force  is  expressed  as 


where: 

ff  =  equivalent  follower  force  load  per  unit  length 

f^ijj  =  electromagnetic  follower  force  (reaction  load  placed  on  gun) 

Xg^  =  axial  position  of  electromagnetic  follower  force 

|>/|  <<  1  for  all  0  <  X  ^  L 

It  can  clearly  be  seen  from  this  form  of  the  follower  force  loading  that  the  distributed 
loading  of  the  barrel  is  proportional  to  the  curvature  at  any  point  along  the  barrel  behind 
(Loading  ahead  of  is  indistinguishable  from  gas  recoil  loads  and  is  accounted  for  in  equation, 
(1),  which  is  based  on  the  rigid  body  acceleration  due  to  all  combined  launch  loads).  When  the 
projectile  is  positively  accelerated,  the  reaction  force  is  negative.  Thus,  at  locations  of  positive 
curvature  (concave  up)  the  equivalent  transverse  loading  is  negative,  exacerbating  the  deflection. 
The  converse  is  also  true;  locations  of  negative  curvature  are  forced  upwards.  It  is  clear  that 
compressive  electromagnetic  loading  further  aggravates  locations  of  greater  curvature,  behind  the 
accelerator,  in  the  barrel. 

TOTAL  TRANSIENT  SOLUTION 

The  solution  method  used  in  this  study  is  called  the  “Uniform  Segments  Method”  (USM), 
which  was  developed  at  Benet  Laboratories  in  the  mid-1980s.^^  It  employs  a  modal  analysis 
technique  in  which  the  gun  tube  is  sectioned  into  a  number  of  uniform  segments  within  which  the 
Euler-Bemoulli  beam  equation  for  free  vibration  is  applied.  (For  elaboration  on  Euler-Bemoulli 
beam  equations  and  modal  analysis  techniques,  see  References  [8],  [9],  or  [10]).  The  boundary 
conditions  of  the  beam  are  free-free,  and  intersegment  continuity  of  displacement,  slope,  bending 
moment,  and  shear  force  is  preserved  at  the  interfaces.  Each  segment  possesses  a  unique  mode 
shape  function  for  each  natural  frequency.  These  functions  satisfy  the  continuity  requirements 
cited  above  and  are  orthogonal.  The  terms  of  these  functions  contain  the  usual  trigonometric  and 
hyperbolic  forms  seen  in  uniform  beam  analysis.  The  natural  frequencies  and  coefficients  of  the 
terms  are  solved  from  the  master  system  matrix  formed  from  the  boundary  and  intersegment 
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continuity  conditions.  Loads  due  to  support  reactions  and  nonstructural  mass  (such  as  the 
breech)  are  superimposed  with  the  load  formulations  of  equations  (1)  through  (4)  and  applied  as 
external  loading  during  the  transient  solution.  They  are  not  accounted  for  in  the  modal  analysis 
portion. 


Using  modal  analysis  techniques,  the  dynamic  profile  of  the  bore  center  line  may  be 
expressed  as  the  product  of  two  single  variable  functions,  one  in  space  and  one  in  time.  This  is 
the  method  of  separation  of  variables  in  differential  equations 

4^)  (5) 


where: 

N  = 

number  of  mode  shapes  in  analysis 

M 

number  of  uniform  segments  in  analysis 

<Pij(x)  = 

mode  shape  function:  j*  segment,  i^’  mode 

qi(t)  = 

Note:  <Pjj(x)  may  only  be  nonzero  for  x  within 
the  segment 

displacement  amplification  coefficient:  i*  mode 

This  method  effects  a  basis  transform  from  the  time  evolution  of  an  infinite  number  of 
transverse  displacements,  corresponding  to  the  continuum  along  x,  to  the  time  evolution 
(amplification  coefficients)  of  an  infinite  number  of  orthogonal  mode  shapes.  Later  in  this  paper, 
justification  is  provided  for  truncating  the  number  of  natural  vibration  mode  shapes  to  N.  This 
provides  the  motivation  for  the  basis  transform.  The  advantage  of  USM  over  other  modal 
techniques  is  its  piece-wise  analytic  formulation  of  the  mode  shapes,  which  lends  itself  to  the 
determination  of  the  curvature-dependent  loadings  of  equations  ( 1 )  through  (4).  Standard 
numerical  procedures  for  modal  analysis  vibration  solutions  are  used  in  solving  for  the 
amplification  coefficients.  For  a  detailed  exposition  of  USM,  see  Reference  [11]. 

DEVELOPMENT  OF  HYBRID  GUN  VIBRATIONAL  MODEL 

The  two  main  concerns  in  developing  a  gun  model  for  use  by  USM  are  that: 

•  The  modal  parameters  of  the  model  accurately  represent  the  modal  parameters  of 
the  component. 

•  A  minimum  number  of  mode  shapes  be  selected  to  establish  model  convergence 
before  any  reliable  calculations  are  made. 

These  concerns  are  addressed  below. 
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Establishment  of  Modal  Parameters 

To  model  transverse  vibrations  of  a  gun  via  the  modal  technique  of  equation  (5)  requires 
an  accurate  representation  of  the  natural  frequencies  and  mode  shapes  of  the  barrel.  This  entails 
the  selection  of  the  boundary  locations,  normalized  weight,  and  bending  resistance  of  each 
segment,  such  that  the  free  vibration  frequencies  and  mode  shapes  of  the  entire  model  accurately 
mimic  the  component.  Mode  shapes  and  frequencies  determined  by  a  finite  element  method 
(FEM)  implemented  in  MATLAB  were  used  to  provide  a  comparative  set  of  modal 
parameters.'*^^ 

The  geometric  and  model  representations  of  the  gun  tube  are  shown  in  Figure  1.  The 
geometry  is  simple,  because  the  system  was  designed  for  experimental  validation.  The 
MATLAB  model  employed  a  total  of  74  beam  elements,  whereas  the  USM  model  has  three 
segments  that  were  chosen  so  that  the  total  beam  mass  was  preserved. 


FEM  Model:  74  Elements 


0.1 

B  0 

-0.1 


Figure  1.  Cylindrical  barrel  geometry  and  spatial  discretization  of  FEM  and  USM 
m^els.  Hash  marks  represent  nodes  between  finite  elements  and  uniform  segments, 
respectively. 
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The  first  six  natural  frequencies  are  shown  in  Figure  2  along  with  vibrational  mode 
shapes.  The  modes  match  well  with  an  increasing  difference  in  the  higher  frequencies.  In 
general,  the  method  of  modal  analysis  simplifies  the  dynamic  solution  because  arguments  based 
on  structural  damping  and  mechanical  impedance  imply  that  at  some  — application  specific — 
cut-off  frequency,  the  inclusion  of  higher  modes  contributes  a  negligible  amount  to  the  dynamic 
solution.  It  has  been  shown  that  eight  mode  shapes  (two  rigid  body,  and  six  flexible)  are 
sufficient  to  model  the  60-mm  gun  barrel.^*^^ 

Mode  Shape  1 ,  FEM  27.48Hz,  USM  26.4  Hz  Mode  Shape  4,  FEM  338.2Hz.  USM  329.5  Hz 


Normalized  Axial  Position  Normalized  Axial  Position 

Figure  2.  Mode  shape  and  frequency  comparison  of  the  barrel  computed  using 
FEM  {dotted  line)  and  USM  {solid  line).  Mode  shapes  have  been  unit-normalized 
in  both  axial  length  and  transverse  deflection. 


RESULTS 

The  gun  dynamic  simulation  of  a  hybrid  launch  was  run  on  the  cannon  designed  and 
constructed  at  Ben6t  Laboratories  in  collaboration  with  an  electrical  engineering  group  from 
Polytechnic  University  at  Brooklyn,  The  simulations  were  conducted  for  an 

armature/projectile  of  mass  412  g,  with  46.5  g  of  IMR4227  propellant.  This  provides  a  gas- 
propelled  speed  of  about  600  m/s.  For  this  simulation,  the  electromagnetic  augmentation  was 
provided  by  a  two-coil  accelerator  configuration  with  the  coils  separated  by  a  distance  equal  to 
the  width  of  two  coils.  (A  full  configuration  for  this  system  consists  of  12  coils  operating  in  3 
phases.  This  two-coil  configuration  provides  for  initial  testing  of  one-half  of  a  circuit  for  a  single 
phase.)  The  two  6-tum  coils  were  powered  by  a  150  pF  capacitor,  with  an  initial  voltage  of  5 
kV.  The  electromagnetic  force  trajectory,  shown  in  Figure  3,  was  computed  using  the 
electromagnetic  finite  element  package,  Flux-2D.‘’^^ 
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Time  (sec)  ^  jq-j 

Figure  3.  Electromagnetic  force  applied  to  the 
armature/projectile.  The  force  on  the  barrel  is  the 
negative  of  this. 


The  dynamic  flexure  of  the  hybrid  launch  is  displayed  in  Figure  4  for  three  different 
mounting  conditions.  As  indicated  in  the  title  of  the  subplots,  the  location  of  the  forward  of  two 
hangers  is  located  at  three  locations;  1.3, 2.5,  and  3.7  meters,  while  the  rear  hanger  is  maintained 
at  0.4  meter.  The  hangers  support  the  gun  system,  effecting  a  simplified  gun  mount  that  allows 
static  deflection  of  the  barrel.  The  plots  depict  the  dynamic  muzzle  pointing  angle.  (Note  that 
this  does  not  include  the  static  contribution).  Inspection  reveals  that  the  deflection  oscillates  as  a 
function  of  time,  and  that  the  modulation  (phase  and  frequency)  of  the  oscillation  may  be  altered 
via  the  parametric  changes  to  the  system. 

The  effect  of  the  electromagnetic  loading  in  these  cases  may  be  seen  as  the  discontinuity 
that  occurs  near  0.007  second.  In  the  first  case,  1.3  meter  hanger  position,  the  loading  improved 
the  performance  of  the  gun  system  by  giving  the  muzzle  a  slight  straightening  impulse. 

(Although  the  loading  of  Figure  3  is  not  truly  impulsive,  it  is  very  quick  relative  to  the  firing 
cycle.)  In  the  second  case,  the  same  effect  is  much  less  pronounced  due  to  the  low  angular 
deflection  of  the  muzzle  (and  the  accelerator).  In  the  third  case,  the  loading  drove  the  muzzle 
further  away  from  its  static  deflection. 

It  is  desirable  to  design  the  system  with  the  least  amount  of  angular  deflection  at  shot 
ejection  to  minimize  round-to-round  dispersion.  Achieving  this  design  objective  is  complicated 
in  practice  by  the  different  interior  ballistics  of  various  rounds.  A  design  that  is  well  timed  for 
one  round  may  be  poorly  timed  for  alternate  rounds.  Thus,  a  compromised  design  is  sought  to 
reduce  variations  and  sensitivity  across  a  wide  range  of  possible  round  types. 

DOWN-BARREL  RECOIL  FORMULATION 

The  purpose  of  this  section  is  to  explicitly  derive  the  transverse  beam  loading  that  results 
from  recoil  loads  that  do  not  necessarily  occur  at  the  zero  axial  position.  Such  loads  may  be 
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Time  {sec} 


Time  (sec)  jq-3 


Figure  4.  Hybrid  gun  dynamic  simulation  for  three  different  system 
parametric  conflgurations.  The  three  main  plots  depict  the  dynamic 
contribution  to  the  angular  deflection  of  the  barrel,  at  an  axial  position 
of  3.7/m,  versus  time.  The  insets  depict  the  dynamic  contribution  to  the 
bore  profile  at  shot  ejection. 
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imposed  by  down-barrel  electromagnetic  accelerators.  This  derivation  will  develop  the  method 
used  by  Simkins,^'*^  use  the  same  assumptions/approximations,  and  extend  it  to  the  current 
application. 


The  Method  of  Simkins.  Equations  (1)  throut 


•  Development  of  Simkins  Equation  1 


“During  the  recoil  of  a  gun  tube  there  results  an  axial  load  per  unit  length  which  is  equal 
to  the  product  of  the  recoil  acceleration,  a(t),  and  the  mass  per  unit  length  of  the  tube, 
p(x).”  This  may  be  interpreted  as  a  distributed  D'Alembert  inertial  force  loading 


=  -p{x)a{t) 


•  Development  of  Simkins  Equation  2 

“When  the  tube  is  curved,  this  load  creates  a  moment  at  any  location  x,  along  the  tube.” 


This  integral  essentially  integrates  the  distributed  D'Alembert  loading  times  the  relative 
lateral  moment  arm  for  all  points  beyond  the  point  of  interest,  x.  One  may  qualitatively 
reason  that  at  a  given  axial  position,  x,  the  effect  of  all  barrel  inertial  beyond  x,  passes 
through  the  barrel  at  x.  Thus  the  integration  from  x  to  L.  (This  formulation  is  in 
agreement  with  the  Extended  Hamilton  Principle  Formulation.)'*^^ 

•  Development  of  Simkins  Equation  3 

Under  the  Euler  beam  assumption  that  all  strain  energy  is  contained  within  the  bending 
mode,  excluding  shear  and  axial  strain,  the  distributed  moment  of  equation  (7)  may  be 
recast  into  an  equivalent  distributed  transverse  loading.  (One  that  would  result  in  the 
same  distributed  moment  as  equation  (7).)  This  distributed  transverse  loading  is  the 
second  spatial  derivative  of  equation  (7)  under  the  small  angle  approximation. 
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The  result  obtained  is 


(8) 


Note  that  this  result  is  equivalent  to  equation  (1)  for  the  uniform  mass  density  within  each 
segment  empolyed  by  the  USM.  (A  slight  unit  ambiguity  exists  in  the  compact  notation 
of  equation  (1)  that  is  better  clarified  by  equation  (8).  The  integral  in  equation  (8)  adds  a 
length  unit  to  the  numerator  of  the  second  term  that  is  missing  in  equation  ( 1 ).) 

•  Explicit  Derivation 

The  double  spatial  differentiation  of  equation  (7)  to  determine  equation  (8)  is  obtained  as 
follows. 

First,  invert  limits  of  integration  of  equation  (7),  insert  equation  (6),  and  separate  the  two 
components  of  the  integral 


X 

M(x,r)=Jw(Ti,r)[y(Ti,r)-y(x,r)]Jri  (9a) 

L 

X 

=  J-p(il)a(t)[y(ii,r) -y{x,t)]dr\  (9b) 

L 


L  L 


Second,  twice  differentiate  equation  (9c)  with  respect  to  x,  recalling  the  chain  rule  and 
that  the  differential  of  a  function  under  integration — with  respect  to  its  upper  limit  of 
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integration — is  the  function  itself,  evaluated  at  the  upper  limit 


X 

M'{x,t)  =  a(r)y(x,r)p(x)  +  a(t)y'(x,t)Jp(ri)dr] 

L 

-  a{t)y[x,t)p{x) 

X 

y  '{x,  t)  p(x)  +  y  "{x,  t)Jp(Tj)dr] 


M%x,t)=a{t) 


M"{x,t)  =  a(/) 


L 

L 


-y'U>f)Jp(n)dr] 


(lOa) 


(10b) 


(10c) 


D'Alembert  Force  Balance 

It  is  important  to  realize  that  the  spatial  integral  of  w(x,t)  along  the  entire  barrel  is  equal 
and  opposite  to  the  net  external  axial  force  placed  upon  the  barrel.  For  gas  launch,  this  load  is 
solely  attributable  to  the  propellant  pressure  at  the  breech.  (Friction  of  the  projectile  and  other 
high-order  contributions  are  generally  neglected.  Also,  many  modeling  techniques,  such  as  in 
Reference  [7],  treat  extraneous  mass  attached  to  the  barrel  as  an  external  loading.)  The 
D'Alembert  force  balance  may  be  written 


(11a) 


(11b) 


where: 


K 

AXIAL  (x,t)  = 

FD(X’t) 


number  of  external  axial  load  sources 
distributed  axial  loading  for  the  d***  load  source 
total  distributed  external  loading 
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Gas  Launch  Contribution 


The  gas  launch  contribution,  neglecting  minor  bore  diameter  variations  within  the 
chamber,  may  be  listed  as 

(12a) 
(12b) 


//(O  = 

/i^(0  =  -■P(0^i,sin|y'(x,r)|^  j6{x-0) 


The  first  term,  the — cosine — component,  is  the  force  that  is  parrallel  to  the  undeformed 
axis.  The  perpendicular — ^sine — component  may  be  directly  included  in  the  transverse  loading. 

Down  Barrel  Axial  Loading 

The  down  barrel  electromagnetic  follower  force  may  be  expressed  in  the  context  of 
equation  (1 1)  as 


/aCO  = /JO  cos 

f/(^,o|  ] 

/2'(0  = /JO  sin 

V(^,ol  ] 

(13a) 

(13b) 


(Again,  the  perpendicular — sine — contribution  must  be  directly  incorporated  into  the  transverse 
loading.) 
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Distributed  Moment  Due  to  Axial  Loading 


In  analogy  with  equations  (7)  and  (10),  the  external  axial  loads  generate  distributed 
moments  along  the  deformed  barrel  expressible  as 

L 

M(x,r)  =|-F^(r|,r)[y(Ti,r)  -y(x,/)]<fTi  (14a) 

X 

L 

M  "{x,  t)=y  '{x,  t)  ( -Fp(x,  t))  -  y  "{x,  t)J(-  FJt},  t)jdr]  (14b) 

Jt 


Distributed  Lateral  Loading  Due  to  External  Axial  Loads 

Equation  (14b)  represents  the  equivalent  distributed  transverse  loading  of  all  combined 
external  axial  loads.  In  both  of  the  loads,  equations  (12)  and  (13),  the  unity  integration  property 
of  the  Dirac  6  function  is  used  to  compute  the  second  term  of  equation  (14b).  Also  note  that  the 
incoporation  of  a  truely  distributed  axial  loading — such  as  would  be  obtained  by  a  long  traveling 
wave  multi-coil  accelerator  stack — is  directly  enabled  by  equation  (14b). 

•  Lateral  Loading  Due  to  Gas  Launch 

The  effective  lateral  loading  due  to  the  pressure  at  the  breech  may  be  expressed  as  the 
combination  of  equation  (12b)  and  the  second  spatial  derivative  of  the  moment  created  by 
equation  (12a)  as  shown  in  equation  (14b) 

/j(x,r)  =  -F(r)A^sin^y'(x,0|^  j6(x-0) 

+  3't^=0(-[-^(0^fcCOs^y'(x,r)|^  J  6(x-0)jj 

(15) 

-y"(x,t)||  -J-F(r)A^cos|y'(ti,r)|^ 


In  this  expression,  the  first  term  is  the  direct  perpendicular  component  of  the  breech 
force,  equation  (12b).  The  second  term  cancels  the  first  under  the  small  angle 
approximation.  The  third  term  only  has  a  non-zero  value  at  x  =  0.  By  the  definition  of 
the  Dirac  function,  the  integral  has  the  value  of  the  function  within  it  at  q  =  x  =  0. 
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However,  the  resultant  force  distribution  is  finite  over  an  infinitesimal  length — a  null 
function;  further  the  boundary  condition  of  the  free  end  of  a  beam  imposes  zero  curvature 
at  X  =  0.  Thus  the  entire  expression  of  equation  (15)  contributes  negligibly  to  the  gun 
dynamics  solution.  (Note:  small  angle  assumptions  are  intermixed  with  the  Euler  beam 
approximations.  Perturbation  and  stability  analysis  would  require  a  detailed  tracking  of 
all  assumptions  to  prevent  confusion  with  order  approximations.) 

Lateral  Loading  Due  To  Electromagnetic  Launch 

The  effective  lateral  loading  due  to  the  electromagnetic  accelerator  may  be  expressed  as 


L 

-y"M J  I  -^Jt)cos|  j6(Ti  -xjjj  dr\ 


(16) 


As  in  the  previous  case,  the  first  two  terms  cancel  under  the  small  angle  approximation. 
However,  the  third  term  effects  a  lateral  loading  as  long  as  x^^  resides  within  the  limits  of 
integration;  this  may  be  seen  to  result  in  the  Heaviside  function.  This,  combined  with  the  sign 
cancellation  and  small  angle  approximation  for  the  cosine  term  results  in  equation  (4). 

The  appearance  of  the  Heaviside  function  may  create  confusion  when  the  electromagnetic 
accelerator  is  decelerating  a  round.  (This  may  be  done  to  tune  the  velocity  of  a  round,  or  may 
occur  during  normal  operation  as  can  be  seen  in  Figure  3  near  7.08  seconds.)  The  confusion 
arises  due  to  the  reversal  of  the  problem;  the  affect  is  to  compress  the  barrel  between  the  axial 
location  of  the  accelerator  and  the  muzzle.  Since  the  Heaviside  is  zero  in  that  region,  it  seems  as 
though  the  expression  does  not  properly  address  this  mode. 

The  confusion  is  eliminated  by  realizing  that  the  electromagnetic  contribution  to  the  gross 
acceleration  of  the  entire  barrel  is  included  within  the  recoil  inertia  load  of  equations  (1)  and  (8). 
Li  fact,  as  can  be  seen  in  equation  (11b),  the  external  axial  loads  and  D'Alembert  inertial  load 
could  be  combined  prior  to  taking  the  second  spatial  derivative  of  the  effected  moment  in 
equations  (8),  (10),  and  (14).  Equation  (4)  may  be  thought  of  as  undoing  the  transverse  inertial 
loading  component  that  is  attributable  to  the  accelerator  force,  between  the  breech  and 
accelerator,  when  in  the  decelerating  mode  of  operation. 
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CONCLUSIONS 


The  method  presented  in  this  paper  will  address  accuracy  considerations  early  in  the 
design  of  electromagnetically  augmented  gun  systems.  Simulation  data  may  be  used  to 
parametrically  optimize  the  design  of  the  gun  system  from  a  structural  dynamics  perspective. 
This  is  particularly  relevant  for  hybrid  gun  systems.  Unlike  gas  propulsion,  whose  recoil  force  is 
manifest  at  the  breech,  pulling  the  inertia  of  the  gun  barrel,  hybrid  electromagnetic  recoil  forces 
are  located  near  the  muzzle,  pushing  the  inertia  of  the  barrel.  This  may  exacerbate  the  affect  of 
gun  whip  but  can  largely  be  avoided  by  altering  the  system  parameters  to  time  the  augmenting 
force  to  coincide  with  low-vibration  amplitudes  of  the  fundamental  modes  of  the  barrel. 
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